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Abstract 

The topographical scattering of gravity waves is investigated using a spectral 
energy balance equation that accounts for first order wave-bottom Bragg scat- 
tering. This model represents the bottom topography and surface waves with 
spectra, and evaluates a Bragg scattering source term that is theoretically valid 
for small bottom and surface slopes and slowly varying spectral properties. The 
robustness of the model is tested for a variety of topographies uniform along one 
horizontal dimension including nearly sinusoidal, linear ramp and step profiles. 
Results are compared with reflections computed using an accurate method that 
applies integral matching along vertical boundaries of a series of steps. For small 
bottom amplitudes, the source term representation yields accurate reflection es- 
timates even for a localized scatterer. This result is proved for small bottom 
amplitudes h relative to the mean water depth H . Wave reflection by small 
amplitude bottom topography thus depends primarily on the bottom elevation 
variance at the Bragg resonance scales, and is insensitive to the detailed shape 
of the bottom profile. Relative errors in the energy reflection coefficient are 
found to be typically 2h/H. 
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Surface waves, scattering, wave reflection, spectral analysis. 
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INTRODUCTION 



Wave propagation over any bottom topography can now be predicted with 
boundary element methods or other accurate numerical techniques. However, 
wave forecasting relies to a large extent on phase-averaged spectral wave models 
based on the energy or action balance equation (Gelci et al. 1957) . For large 
bottom slopes waves can be reflected and this reflection is currently not repre- 
sented in these models, while the significance of this process is still poorly known 
(Long 1973; Richter et al. 1976; Ardhuin et al. 2003). For waves propagating 
over a sinusoidal seabed profile, a maximum reflection or resonance is observed 
when the seabed wavenumber is twice as large as the surface wave wavenumber 
(Heathershaw 1982). Davies and Heathershaw (1984) proposed a deterministic 
wave amplitude evolution equation for normally incident waves over a sinusoidal 
seabed, based on a perturbation expansion for small bottom undulations. This 
theory was shown to be in good agreement with experimental data but overes- 
timates reflection at resonance. Mei (1985) developed a more accurate approx- 
imation that is valid at resonance using a multiple scale theory. This approach 
was further extended to random bottom topography in one dimension (Mei and 
Hancock, 2003). The Bragg resonance theory can be extended to any arbi- 
trary topography in two dimensions, that is statistically uniform(Hasselmann 
1966). Ardhuin and Herbers (2002) further included slow depth variations. 
The resulting spectral energy balance equation contains a bottom scattering 
source term Sbscati which is formally valid for small surface and bottom slopes 
and slowly varying spectral properties. Sbscat is readily introduced into exist- 
ing energy-balance-based spectral wave models, and was numerically validated 
with field observations (Ardhuin et al. 2003). While this stochastic theory is 
in a good agreement with deterministic results for small amplitude sinusoidal 
topography (Ardhuin and Herbers 2002), the assumed slowly varying bottom 
spectrum is not compatible with isolated bottom features, and the limitations 
and robustness of the source term approximation for realistic continental shelf 
topography arc not well understood. The limitations of the stochastic source 
term model are examined here through comparisons with a deterministic model 
for arbitrary one-dimensional (ID) seabed topography that is uniform along the 
second horizontal dimension. We review the random Bragg scattering model, 
and investigate the applicability limits of the source term for a variety of seabed 
topography. Predicted reflection coefficients are compared with results based on 
Rey's (1992) model, which approximates the bottom profile as a series of steps. 
Examples include modulated sinusoidal topography that is well within the va- 
lidity constraints of the source term approximation as well as a steep ramp and 
a step that violate the assumption of a slowly varying bottom spectrum and 
thus provide a simple test of the robustness of the source term approximation. 
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THEORETICAL BACKGROUND 
MATCHING BOUNDARY SOLUTION 

We use Rey's (1992) algorithm, based on the theory of Takano (1960) and Kirby 
and Dalrymple (1982). It uses a decomposition of the bottom profile in a series 
of ./V steps with integral matching along vertical boundaries between each pair of 
adjacent steps. A coordinate frame is defined with the horizontal x coordinate 
in the direction of the incident waves and the vertical z coordinate pointing 
upwards relative to the mean water level. The velocity potential is described 
by a sum of flat bottom propagating and evanescent modes. Evanescent modes 
are included in the matching condition to ensure a consistent treatment of the 
wave field (Rey 1992). The general solution of the velocity potential for a step 
(p) of depth Hp is given by the following equations: 

%{x,z,t) = 4>p{x,z)e-' lwt for p=l,N, (1) 

with, 

Q 

<t> P {x,z) = A±e ±ikx Xp (z) +J2 B PA e±k * X ^«(^ ( 2 ) 

q=l 



propagating modes 



evanescent modes 



where (x P , ip P q q—l.Q) define a complete orthogonal set for each step region 
(P): 

X P (z) = cosh k p (H p + z), (3) 
ip P , q (z) = COS kp. q (Hp + z). (4) 
kp and k pq satisfy the following dispersion relations: 

up" 

— — k p tanh(fcpiip), (5) 

-j- = -fcp, 9 tan(fc p , 9 i/ p ). (6) 
where g is the acceleration of gravity. 

Across each step (p), matching conditions between two domains (labelled 
p = 1 and p = 2 in figure ^) must be applied to ensure continuity of the fluid 
velocity and surface elevation. 

* = 7£ = 7j£f fOT ~H 2 <z<0, (7) 

^fi = for -H Y <z< -H 2 . (8) 
ox 
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Figure 1: Stepwise approximation 



The integral formulation of these conditions (for Hi > H 2 ) leads to: 
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(9) 
(10) 

(11) 
(12) 



The orthogonality of the set functions largely simplifies these equations. In 
order to solve the problem numerically, the number of evanescent modes q are 
truncated to q — Q. Practically, only a few evanescent modes are needed to 
ensure convergence. For N steps, 2N(Q + 1) equations are solved to obtain 
the 2N(Q + 1) complex coefficients and . At the boundaries (p=0 and 
p=N), the reflection coefficient is given by: 



(13) 



This method has the advantage that it is valid for arbitrary ID topography. 



BRAGG SCATTERING THEORY 

We consider random waves propagating over a 2D irregular bottom with a slowly 
varying mean depth H and small-scale topography h. The bottom elevation is 
given by z = —H(x.) + h(x), with x the horizontal position vector. The free 
surface position is ((x,t). 
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Figure 2: Definitions 



Considering an irrotational flow for an incompressible fluid, we have the gov- 
erning equations and boundary conditions for the velocity potential <fi: 

V 2 </>+^$ = for -H + h<z<(, (14) 
oz z 

^ = V.4> ■ Y(/i - H ) at z = -H + h, (15) 
oz 

m = Tz at z = c (16) 

9C+f t =-\[\Y_4>\ 2 + ( d -^) 2 ] at Z = C, (17) 

where V and V 2 are the horizontal gradient and Laplacian operators. The 
equations <|14|l . I)15[l. Hlfijl and Ijl7|l are respectively the Laplace's equation, free 
surface and bottom boundary conditions, and Bernoulli's equation. Combining 
these two last equations, we obtain: 

d 2 4> d4> dW<j> d<j> d 2 4> . , 

Assuming that the surface and the small-scale bottom slopes are of the same 
order e, and the large scale bottom slope is of order e 2 , a perturbation expansion 
of <fi up to the third order in e yields the following spectral energy balance 
equation (details are given in Ardhuin and Herbers 2002) : 

dE(k,x,t) _ 

3£ — <Jbscat(k,x,*)j I. -!■»,) 

where 

/■27T 

S bscat (k, x,t)=K{k,H) / cos 2 (0-0') (k - x) [-E (k^, x.,t) — E (k, x, i)] d0', 

Jo 

(20) 

with 

47T[jfc 4 

^ ^' ^ = sinh (2kH) [2kH + sinh (2fcff )] ' ^ 
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E(k,x,t) is the surface elevation spectrum and F B (k, x) is the small-scale bot- 
tom elevation spectrum. These spectra are slowly varying functions of (x, t) 
and x resectively. k is the wavenumber vector defined by k = (k cos 9, k sin 9) = 
(k x , k y ), where 9 defines the angle with the x-axis. The spectral densities E and 
F B are defined such that the integral over the entire k-plane equals the local 
variance, 

/ + OO I- + CC 
/ F B {k,x)dk x dk y . (22) 

The frequency to is given by the dispersion relation: 

u? = gktaii\i{kH). (23) 

Here we consider a steady wave field in one dimension with incident and reflected 
waves propagating along the x-axis. After integration over k y , k x becomes k 
and (|T$)l reduces to 



with a source term 



F B (2k t) 

S hscat (k 1 x)=K(h,H) L-U.[E(-k,x)-E{k,x)]. (25) 



The first term of Eq. (|24|l represents advection in physical space with the group 
velocity defined by 

c >= d i=t <-) 

and the second term describes the effect of shoaling on the wavenumber 



= dk 2fc 2 dH 

k dt 9 2kh + sinh(2fc/i) dx 



REFLECTION BY MODULATED SINUSOIDAL 
BOTTOM TOPOGRAPHY 

The source term approximation was validated by Ardhuin and Herbers (2002) 
for random waves reflecting from a sinusoidal seabed, by integrating 5bscat an- 
alytically across the wave spectrum in the limit of weak reflection (E(—k) << 
E{k) , with positive and negative wavenumbers corresponding to the incident 
and reflected waves, respectively). A comparison with Dalrymple and Kirby's 
(1986) solution gave good agreement, even for only a few bars. For stronger 
reflection, equation l|24|l is not readily evaluated analytically, and numerical 
integration is not feasible since a sinusoidal bottom has an infinitely narrow 
spectrum (a Dirac distribution), and thus cannot be represented with a finite 
bottom discretization Akf,. 



6 



0.42 
0.415 

0.41 
0.405 

E 

I M 
D 

0.395 
0.39 
0.385 



""-5 -4 -3 -2 -1 1 2 3 4 5 

x(m) 

Figure 3: Modulated seabed (m=3), 6fc& o = 0.06 

We consider instead a bottom spectrum with a finite width that corresponds 
to a modulated sinusoidal bottom profile. 

The modulated seabed is represented by a sum of cosines: 

h(x) = S^~i ) 1 / j , /2 6 i aw[(A5,o + iAk b )x] (28) 

The slowly varying depth (H ) , defined in part 2 is taken constant while the per- 
turbation (h) represents the modulated seabed. We define the root mean square 
(r.m.s.) bar amplitude b from the bottom variance, b = v 7 < h 2 >, and a repre- 
sentative bottom slope e = bkb,o- The reflected wave energy is calculated for the 
bed profile shown in figure with the peak bottom wave number fc^o = 6m _1 
(•V>.o = 1.04m), and a short modulation length with m = 3, and equal ampli- 
tudes (bi) for all bottom components. The length of the bed is 1.5 modulation 
lengths, giving the bottom spectrum shown in figure The reflection from this 
modulated sinusoidal bottom was evaluated for an incident Pierson-Moskowitz 
spectrum, with a peak at fco satisfying the Bragg resonance condition 2fco = &/>.o 
(FigEJ). Spectral results for Rey's model were obtained by evaluating reflection 
coefficients for monochromatic waves over a range of frequencies and integrat- 
ing the reflected energy across the spectrum. 70 steps are used to resolved the 
bathymetry. Results for various values of b are displayed in the form of reflection 
coefficients R (FigEJ as a function of the slope bkb,o- R is defined by the ratio of 
the reflected and incident energies: R= (E/ C<0 i?)/ (£fc>o-E)- Predictions based 
on the source term method (i?smod) and the matching boundary model using 
5 evanescent modes (i?MBmod) agree well over a wide range of bottom slopes. 
The solutions gradually diverge for large bottom slopes where the source term 
underpredicts the reflection. Even for bko = 0.3 (6/Aq = 0.05), differences are 
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Figure 4: Modulated seabed spectrum (m=3) 

less than 10% confirming the robustness of the source term method for steep 
topography. 

To evaluate the effect of the spectral width on the reflection coefficient, fig- 
ure also includes predictions for sinusoidal topography (m=0) with the same 
variance. Results for sinusoidal topography were obtained using Mei's (1985) 
analytical approximation and Rey's (1992) algorithm. The resulting reflection 
coefficients i?Moi and i?MBsin, respectively, agree for small bottom slopes (FigJfjJ 
and diverge for larger slopes as already shown by Rey (1992). Indeed, i?Moi was 
derived for small bottom slopes while the matched boundary solution converges 
to the exact reflection for any bottom profile when the number of evanescent 
modes goes to infinity. What may seem surprising is that the reflection coef- 
ficient for the sinusoidal and modulated sinusoidal topographies i?MBmod and 
-RiviBsin agree for small slopes although bottom profiles are quite different. Ap- 
parently, for small bottom slopes and narrow bottom spectra the reflection is 
only a function of the total bottom elevation variance b 2 and does not depend 
on the phases of its components. This result is obvious from the viewpoint of 
the source term theory that was derived for small bottom slopes, and does not 
retain the phases of the bottom spectrum components. The predicted reflection 
depends on the convolution of the wave spectrum with the bottom spectrum at 
the Bragg resonance wavenumber (the integral of 112511 over all wavenumbers) . 
If the bottom spectrum is narrow compared with the wave spectrum then the 
total source term depends only on the total bottom variance and the surface 
spectral density at the Bragg resonance wavenumber. 
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Figure 5: Incident wave spectrum 

REFLECTION BY A LINEAR RAMP 

To investigate the robustness of the variance-based source term model for re- 
flection induced by localized topography, we consider the linear ramp problem 
used in previous studies to test the mild slope equation (Booij 1983). In the 
source term approximation, wave scattering is the result of interactions between 
surface waves and bottom variations at the scale of the surface wavelength. The 
scattering model is thus based on a decomposition of the topography into a 
slowly varying depth H and a perturbation h (small scale topography), which 
corresponds to a separation between refraction and shoaling that occurs over 
the slowly varying depth H and scattering at these short scales. For practical 
applications, it is desirable to have a perturbation h that is zero outside of a 
finite region, so that the spectrum of h is well defined. Once the two criteria 
that the slope of H does not exceed a given threshold and h is zero outside of a 
region of radius nL are satisfied, the choice of the depth decomposition in h and 
H is fairly arbitrary and does not affect the following results. For simplicity we 
take a piecewise linear function for H(x), so that the perturbation h(x) takes 
the form of a triangular wave f FiglTjl . 

The ramp profile is defined by the fixed water depths Hi, H2, while the ramp 
slope a is varied by adjusting its length 2L (FigQ. To ensure that H{x) is slowly 
varying, 7 has to be small. This is achieved by extending the domain to a length 
2nL with n > 1 (FigEJ). The slope of H is then given by tan 7 = (tan a)/n, 
with several values of n tested below. 
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Figure 6: Wave reflection by modulated sinusoidal bottom 
First test case: small depth change 

We first consider a ramp with a small depth transition from H\ = 0.5m to 
H 2 = 0.3m. The incident wave spectrum is represented by the same Pierson- 
Moskowitz spectrum that was used in the previous section with the peak wavenum- 
ber in deep water k = 3m _1 (FigEJ, so that k Hi = 1.5 and k H 2 — 0.9. In 
order to investigate the source term applicability limits, the linear ramp slope 
tan a is varied from 0.01 to 2.9. For each value of a, several values of 7 are tested, 
with n varying from 5 to 50. The reflection coefficient i?s (source term reflection 
due to the residual) is compared with the " exact" computation Rmb (matching 
boundary algorithm) in figure[B]and the relative error (Rs — Rmb) /Rmb is shown 
in figure In our calculations, for slopes of H such as tana < 0.4, i?s,n=5 is 
within 30% of the exact value Rmb- For larger values of tana, Rs.n=5 decreases 
and tends to zero (Fig[SJ) , while the exact solution Rmb converges to the re- 
flection over a vertical step as tana goes to infinity. The value tana = 0.4 
corresponds to tan7 (= tana/5) equal to 0.08. For larger n the slope of H is 
reduced and Rs, n is valid for a wider range of ramp slopes. 

We notice that for all values of n shown in figure^ the model gives reasonable 
results for tan 7 (= (tan a)/n) up to about 0.08. The ramp slope does not 
appear to be a limiting factor (as it was assumed in the theory). For tan 7 
larger than 0.08 the reflection is increasingly underestimated probably because 
of the contribution of the large scale profile H (x) to the reflection. As n increases 
h approaches the slope of the actual ramp and i?s,n converges to i?s i0 o which 
is about 10% larger than Rmb for all ramp slopes. As discussed below, the 
accuracy of the model is apparently not limited by the ramp slope. 
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Figure 7: Decomposition of linear ramp (solid line) into a slowly varying depth 
H (dashed line) and residual h (dotted line), (a) and (b) for small and large n 
respectively 



It may seem surprising that Rs. n actually converges for large n while the 
bottom spectrum does not. In the case of a vertical step of height h in the 
middle of a domain of length 2nL, the spectral density F B (k) of a discrete 
variance spectrum of the residual is proportional to h 2 /2nLk 2 and tends to 
zero (except around k = 0) as n goes to infinity. However the source term 
formulation represents scattering as uniformly distributed along the bottom, 
and the integration of the source term along the wave propagation path yields 
a reflection that is proportional to 2nL F B (k) and thus converges when n goes 
to infinity. The use of infinite support for H and h (taking the limit n — ► oo) to 
compute the reflection over a localized ramp is counterintuitive. It represents 
a physically localized scattering with a mathematically distributed source. In 
practice, the bottom spectrum is obtained by discrete Fourier Transform of the 
bottom, and it only tends to continuous power spectrum in the limit n — > oo. 
Further, it should be realized that the bottom power spectrum is the Fourier 
transform of the bottom autocorellation function used by Mei and Hancock 
(2003, see Appendix). 

For a non-random bottom such as the ramp here, one may use intermediate 
results by Mei and Hancock (2003) where the hypothesis that the bottom is 
random only comes in for discarding nonlinear wave effects (which are not taken 
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Figure 8: Wave reflection by a ramp 

into account here). It thus appears that our rather surprising result for the 
convergence as n — > oo is justified by the convergence of the discrete spectrum 
to the continuous power spectrum and the theory of Mei and Hancock (2003) 
applied to non-random bottoms (see Appendix). It shows that the far field 
scattered energy by small amplitude depth variations only depends on the power 
spectrum of the scatterers at the Bragg scale, and not on its localization in space, 
as long as the bottom amplitude remains small. 

Booij's ramp: larger depth change 

This approach should clearly break down for finite bottom amplitudes, in partic- 
ular because sub-harmonic scattering was observed (Belzons et al. 1991) while 
it is not explained by the present theory. Such a limit should be tested to see 
whether our present approach has some practical applicability. We therefore 
take a second test case with a larger ramp is taken from Booij (1983) with 
water depths Hi = 4.97m, H2 = 14.92m and an incident wave peak period 
T = 10s. The corresponding peak wavenumbcr in deep water ko — 0.04m _1 so 
that koHi = 0.6 and fco#2 = 0.2. Results for ramp slopes tana ranging from 
0.001 to 2.9, and n=10 and 50 are shown in figure ITUl 

We notice again that Rs,n converges for large n, provided that (tana)/n < 
0.08. However, in this case the relative error is larger than in the first test case, 
up to about 30%. The two tests have the same ramp slopes but different ratio of 
water depths at the edges of the ramp H1/H2 = 3 here versus H1/H2 = 1.7 in 
the previous case. The two cases suggest that the source term is more sensitive 
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Figure 9: Relative errors in wave reflection by a ramp 

to the amplitude than the slope of the bottom perturbation h. Formally, the 
bottom amplitude only appears in the bottom boundary condition l|15Jl . which 
is linearized at z = — H using the following Taylor series expansion; 

d$ h 2 <9 2 $ 

$| z= _ H+/l = §| a= _ H + h^\ z= - H + Y^*=- a + ° {k) - (29) 

Ardhuin and Herbcrs(2002) use a representative length scale l/fco to non-dimensionalizc 
(EH) as, 

<9$ n 2 <9 2 $ 

where z = k$z, r\ = koh, r\ corresponding to the scales that cause wave scat- 
tering. The validity of the Taylor expansion requires that r\ is small and also 
that the first and second derivative of <f> with respect to z are of order 1. In 
this approximation (|30J) is limited by the small-scale slope koh. However one 
may also take Ho as the representative length which leads to the same equation 
with T] = h/Ho, limited then by the water depth ratio h/Ho- The choice 
of the representative length was arbitrary and can be justified only a posteriori, 
by evaluating the scale of variation of $ and thus the magnitude of 9$ jdz and 
d 2 $/dz 2 . The numerical results presented here show that the source term is 
more sensitive to the water depth change h/H than the small-scale slope koh. 
Booij (1983) had found that the standard mild slope equation (Bcrkhoff 1972) 
gave errors less than 10% for tan a up to 1 /3. Our results suggest that the Bragg 
scattering model can be as accurate as the mild slope equation for computing 
reflection, but only for Ah/Ho less than 0.2. 
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Figure 10: Wave reflection by Booij's ramp 



REFLECTION BY A STEP 



Now that the effect oih/Ho is well established, one may question the importance 
of other parameters. We thus evaluate source term predictions of broad and 
narrow surface wave spectra over steps of varying height to gain further insight 
into the limitations of the source term approximation for localized topography. 
Reflection of waves by a rectangular step has been investigated analytically 
and experimentally in numerous studies (Neuman 1965a,b; Miles 1967; Mei 
and Black 1969; Mei 1983 and Rey, Belzons and Guazzelli 1992) and is well 
understood. The step is defined in figure ITT1 where 2L is the step-length, h the 
height and 2nL the size of the entire computational domain. 




Free surface 
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-nL 



-L 
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nL 



Figure 11: Sketch of the step 



14 



Numerical set-up 

The spectral density of the bottom F B (k) is proportional to h 2 /2nLk 2 . Hence, 
integration of the source term along the wave propagation path yields a reflection 
that is proportional to 2nLF B (k) ~ h 2 /k 2 , independent of n. Although the 
domain length has no effect on real waves in the absence of bottom friction, it 
influences the discretization of the bottom spectrum (AA: = 2Tr/2nL), and thus 
it may have an impact on the numerical results. However 2nF B (n) converges 
as n goes to infinity ('Fig ll2|l . so that the domain length does not change the 
results for large enough values of n. A large domain with n — 8 was used here. 



■ - Fx n , n=2 
F B x n , n=6 
F B x n , n= 
E, wide 

- - E, narrow 




Figure 12: Wide and narrow surface wave spectra superposed on the bottom 
spectrum for domain sizes n = 2,6 and 8. The bottom spectrum is rescaled by 
the surface wavenumber (F B {k/2)) to show the resonant bottom and surface 
components. 



The step width (2L) is taken to be half the wavelength of the surface waves 
for a spectrum peak ko p = 0.04m _1 (Lq = 157m) in a water depth of 15m. Two 
different wave spectra are used here (bold lines in figure 112")) : a wide spectrum 
(solid) with a classic Pierson-Moskowitz shape, typical of wind seas, and the 
narrow swell-like spectrum (dashed) with a Gaussian shape. Once the shape 
of wave spectrum is chosen, the solution is a function of three non-dimensional 
variables: the step height h/H, the water depth ko p H, and the relative step 
width kopL. 



Influence of the height of the step 

The accuracy of the source term for a range of non-dimensional step heights h/H 
is evaluated in intermediate and shallow water through comparison with the 
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"exact" matching boundary algorithm fFig ll3[) . Energy reflection coefficients 
are compared for two different water depths, ko p H = 0.1 and ko p H = 0.6, 
representative of shallow and intermediate depths. The incident wave spectrum 
has a Pierson-Moskowitz shape. 




h/H 

Figure 13: Reflected energy computed with the source term (dash-line) and with 
the matching boundary algorithm (full- line), for intermediate depth (ko P H = 
0.6) and shallow water (ko p H = 0.1), and relative error of the source term. 

As expected from previous calculations, the error in the source term increases 
with the step amplitude h/H. For h/H < 0.05 the error in the predicted re- 
flection coefficients is less than 10%. These results provide further confirmation 
that the height of the localized scatterer is a limiting factor for the source term 
computation, but not its slope, which is infinite here, and this result holds for 
very shallow water. 

Influence of the width of the step and the wave spectrum 

Here we consider the dependence of the reflection coefficient on the width of the 
step and the width of the wave spectrum for a small amplitude step {h/H = 
0.02) in shallow water (ko p H = 0.1). The non-dimensional step width ko p L is 
varied, effectively changing the position of the wave spectrum peak relative to 
the bottom spectral peaks (see Fig ll2[) . Results are shown in figures ITU and IT51 
for wide and narrow wave spectra, respectively. 

The same computation is done for the narrow spectrum (Fig Et - 
For both wide and narrow surface waves spectra, the source term yields ac- 
curate results, and the errors do not appear to be sensitive to the width of the 
step. Oscillations in the reflection coefficient with varying k p oL represent an 
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Figure 14: Reflected energy computed with the source term (dash-line) and 
with the matching boundary algorithm (full-line) for a wide wave spectrum. 
The bottom spectrum (F B ) is also indicated, scaled by the normalized resonant 
surface wavenumber to indicate the resonant response (bold dash-line). Other 
parameters are h/H = 0.02 and ko p H = 0.1. 

interference phenomenon that has been described in numerous previous studies. 
When a monochromatic incident wave runs up the leading edge of the step at 
x = —L, it is partly reflected and partly transmitted. As the transmitted com- 
ponent passes the rear edge of the step at x = L, it is again partially reflected 
and partially transmitted. If the reflected waves originating from the front and 
rear edges of the step are in phase we have a constructive interference which 
amplifies the reflection. Conversely, destructive interference occurs if the two 
reflected wave trains are 180 degrees out of phase and cancel out, yielding zero 
reflection. For long waves, maximum reflection occurs when sin 2 2k p oL — 1 (Mei, 
1983), where k p o is the incident wave wave number. This condition is met when: 

2k 0p L = (2n-l)^, n = 1,2,3,-.. (31) 

The corresponding values of k p oL are k p oL = 7r/4 ~ 0.78 , 37r/4 ~ 2.35, 57r/4 ~ 
3.93 • • • These values match with the reflection peaks observed in the figures 
H14I15JI both for the source term and the matching boundary algorithm. In 
the wide spectrum case (Fig ll4f) these oscillations are suppressed and for high 
values of k p oL, the reflection tends to a constant value. Using Bragg scattering, 
this is explained by the fact that in the limit of large step width k p oL the 
wave spectrum is wider than the side lobes of the bottom spectrum (see Fig ll2|l 
and the effects of constructive and destructive interferences for different spectral 
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Figure 15: Same as Figll4l but for a narrow wave spectrum. 



component average out. The reflection coefficient is a convolution of the bottom 
spectrum and the surface wave spectrum, and thus the reflection is insensitive 
to bottom spectral details with scales finer than the wave spectrum width. 

CONCLUSIONS 

Predictions of the scattering of surface waves by bottom topography based on a 
spectral energy balance equation that includes a wave-bottom Bragg scattering 
source term (Ardhuin and Herbers 2002) are compared with exact results based 
on a matching boundary algorithm (Rey 1992). The source term yields accu- 
rate reflection predictions for modulated sinusoidal topography. In the limit 
of small bottom amplitudes h compared to the water depth H , the two models 
yield identical results, confirming that the far-field scattered wave is determined 
entirely by the variance spectrum of the bottom and does not depend on the 
phases of its components. This finding also holds for localized topography, a 
result that can be justified by the approach of Mei and Hancock (2003) using 
their intermediate results for non-random bottoms. In that case, the bottom 
spectrum must be carefully calculated over a large enough domain in order to 
resolve the important bottom scales. Using Discrete Fourier Transforms, one 
may use an artificial gently sloping extension of the area covered by scatterers. 
However, it is found that it also holds for very steep topography, such as a single 
step, for a variety of water depths and wave spectrum shapes, as long as h < H 
is small. In our calculations, relative errors in the energy reflection coefficient 
are found to be typically 2h/H, or h/H for the amplitude reflection coefficient. 
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These results show that the Bragg scattering source term is a reasonably accu- 
rate method for representing wave reflection in spectral wave models, for a wide 
range of small amplitude bottom topographies found on continental shelves. 
The source term approach is also very efficient compared to the elliptic models 
such as proposed by Athanassoulis and Belibassakis (1999). An extension of 
the source term to higher order (e.g. following Liu and Yue 1998) may reduce 
errors for larger values of h/ H, that are shown here to be the limiting factor in 
practical applications. Results for ID bottom profiles are expected to hold for 
practical 2D applications of the source term approximation. 
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APPENDIX. RECONCILIATION OF RANDOM 
AND DETERMINISTIC WAVE THEORIES 

Mei and Hancock (2003) considered the same problem of a wave train propa- 
gating over an arbitrary topography of small amplitude h. In their scaling h 
is small compared to the wavelength 2ir/k 1 but, as discussed in this paper, the 
scaling for the bottom perturbation could also be the mean water depth H. 
These authors further assume that h is a random function that is stationary 
with respect to the fast coordinate x, and introduce a slow coordinate x\ for 
variations in the statistics of h. This two-scale approach is similar to that used 
by Ardhuin and Herbers (2002). Mei and Hancock (2003) obtained an ampli- 
tude evolution equation in which the topography acts as a linear damping with 
a coefficient Pi given by their equation (B8) as 

* = ^(fcq) 2 fc(7(2fc)+7(0)) m 
4 cosh 2 kH(cu 2 H/ g + sinh 2 fcif)' 

where a 2 (xi)j is the auto-correlation function of the bottom topography, decom- 
posed in a slowly- varying local variance a 2 (x\ ) and a normalized auto-covariance 
7. 7 is the Fourier Transform of 7. Although Mei and Hancock's (2003) result 
does not conserve energy (which requires the introduction of higher order terms, 
see Ardhuin and Herbers 2002), it is rather general as far as the bottom is con- 
cerned. The essential difference with Ardhuin and Herbers (2002) is that there 
is no need for a large number of bottom undulations to obtain an expression 
for the scattering, and the "number of undulations" is properly defined by the 
scale over which the auto-covariance goes to zero. 
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Naturally the two theories are consistent, and we can obtain from /3j the 
damping coefficient (3e for the energy, which is twice that for the wave amplitude 
A since d(AA*)/dt = = -(3 E AA*, with A* the complex conjugate of 

A. Re-writing i|32|) one has, 

2fcW(7(2fc)+7(0)) m) 
/B sinh 2fci? [2fci? + sinh 2fcif] ' 1 J 

For a zero-mean stationary process the Fourier transform of the auto-covariance 
function is simply 2ir times the power spectral density F B (e. g. Priestley 1981, 
theorem 4.8.1 p 211), so that, for F B (0) = 0, we get 

^ E sinh2kH [2kH + sinh2kHY v ; 

which is the linear part of the bottom scattering source term <|25|) in one dimen- 



sion, 



S bsca , t (k) = f3 E (E(-k)-E(k)), (35) 



Interestingly the hypothesis of randomness for h is not important for the 
value of Pi when averaged over the entire field of scatterers (however, it does 
impact the real part, i.e. the phase of the waves). Following Mei and Hancock's 
(2003) derivation, one may define a that is also a function of the fast coordi- 
nate x using their equation (2.36), and in that case the derivation is identical, 
replacing a 2 (xi)j by h(x)h(x — £), all the way to their equations (B1)-(B3). 
Then one may define a mean value, which, in the case of a finite region with 
scatterers between —nL and nL reads, 

7k = Try [ (3i(x)dx. (36) 
Taking the imaginary part of their equations (B1)-(B3) we have, 

J = i . fc 2 [ Iq ^ kin 1 

1 W 2 cosh 2 (kH) \ uj 2 H/g + sinh 2 (fc J ff) ^ k n [uu 2 H/g + sin 2 (k n H)] J 

(37) 

with (correcting a few minor type-setting errors in their paper), 

I = _R j_L J +OC (JL _ ikj ( h (x)h{x - 0) <- i/£ ' - <lOl r| (38) 

and 



In = -9 ^ II /J - *) {h(x)h(x - 0) e^l^d* 



(39) 
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Switching the order of the integrals. 1381391 are identical to their equations (B2)- 
(B3), provided that we redefine 7 as the full auto-covariance function 



In this case 7 is obviously real and even and we obtain their equation (B8) for 



We have thus proved that in one dimension and in the limit of small bottom 
amplitudes the scattering source term applies to non-random bottoms. In these 
conditions, the linear part of the source term represents the damping of the 
incident waves (and thus also the average scattered wave energy) averaged over 
the area covered by scatterers. 
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APPENDIX.B SYMBOLS 



A = propagating modes amplitude; 

b = root mean square amplitude from the bottom variance; 

B = evanescent modes amplitude; 

Cg = group velocity; 

Ck = spectral advection velocity; 

E = surface elevation spectral density; 

F B = small-scale bottom elevation spectrum; 

h = bottom perturbation height; 

H = water depth; 

k = surface wavenumber; 

fco = peak wavenumber in deep water; 

ko p = peak wavenumber; 

kb,o — peak bottom wavenumber; 

K = Source term coefficient; 

Kr = amplitude reflection coefficient; 

L = half-length of the ramp; 

L = peak wavelength; 

n = mild slope inclination parameter; 

m = modulation parameter; 

R = energy reflection coefficient; 

Rmb = Matching Boundary energy reflection coefficient; 

Ruei = Mei energy reflection coefficient; 

Rs = Source term energy reflection coefficient; 

Sscat — bottom scattering source term for the wave energy spectrum; 

To = peak period; 

a = ramp inclination; 

7 = mild slope inclination; 

s = representative bottom slope; 

( = free surface position; 

r] = small parameter; 

$ = velocity potential; 

X, V'n = complete orthogonal set of functions; 

u> = wave radian frequency; 

Akb = discretization of the bottom spectrum; 

Subscripts 

~= non-dimensionalized variable; 
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